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Abstract 

We present a detailed study of the recently proposed gapless color superconduc- 
tivity in neutral two-flavor quark matter in /^-equilibrium at zero as well as at finite 
temperature. We clarify the nature of the fermionic quasiparticles, and discuss the 
physical properties of this gapless phase of quark matter. The appearance of two 
"secondary" spin-1 condensates in the gapless ground state is pointed out. The gluon 
screening properties are also discussed. It is found that finite temperature proper- 
ties of the gapless color superconductor are rather unusual. One of the most striking 
results is that the ratio of the critical temperature to the gap at zero temperature is 
not a universal number. Moreover, this ratio could become arbitrarily large at some 
values of the coupling constant. In addition, the value of the gap parameter is not a 
monotonic function of temperature, although the phase transition is a second order 
phase transition. In an extreme case, in fact, a nonzero finite temperature gap can 
appear even if its zero temperature value is vanishing. 
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1 Introduction 



Sufficiently cold and dense quark matter is a color superconductor [1]. The 
corresponding ground state is characterized by a condensate of Cooper pairs 
made of quarks. Since the latter carry color charge, the SU(3)c color gauge 
group of strong interactions is partially or completely broken through the 
Anderson-Higgs mechanism. At asymptotic densities, this phenomenon was 
studied from first principles in Refs. [2,3,4]. These studies of dense QCD, as 
well as studies of various phenomenologically motivated models [5,6,7,8,9], 
suggest that dense quark matter at high baryon density has a very rich phase 
structure, consisting of many different phases. 

It is natural to expect that some color superconducting phases may exist in 
the interior of compact stars. The estimated central densities of such stars 
might be as large as lOpo (where po ~ 0.15 fm~^ is the saturation density), 
while their temperatures are in the range of tens of keV. This could provide 
ideal conditions for the diquark Cooper pairing of color superconductivity. 

Matter in the bulk of compact stars should be neutral with respect to electrical 
as well as color charges. Also, such matter should remain in /9-equilibrium. Sat- 
isfying these requirements may impose nontrivial relations between the chem- 
ical potentials of different quarks. In turn, such relations could substantially 
influence the pairing dynamics between quarks, for instance, by suppressing 
some color superconducting phases and by favoring others. For example, it 
was argued in Ref. [10], that the mixture of the 2-fiavor color superconducting 
(2SC) phase and normal strange quarks is less favorable than the color-flavor- 
locked (CFL) phase if the charge neutrality condition is enforced (note that 
the strange quark mass was chosen too small to allow the appearance of a 
pure 2SC phase). A similar conclusion was also reached in Ref. [11]. 

Recently, it was proposed that neutral two- flavor quark matter in /9-equilibrium 
can have a rather unusual, stable ground state called a gapless color supercon- 
ductor (g2SC) [12]. A similar phase of gapless color superconductivity could 
also appear in strange quark matter when the baryon number density is kept 
fixed [13] [for earlier studies see Ref. [14] where an unstable gapless phase 
was considered]. This gapless strange quark matter may, in principle, be real- 
ized inside the so-called strangelets. It was also suggested that closely related 
phases may exist in two component mixtures of cold fermionic atoms [15]. 
It is fair to mention that, in both quark and cold atom systems, the ground 
state could also be given by some mixed phases [16,17,18] provided the surface 
tension effects are not too strong. 

While the symmetry of the g2SC ground state is the same as that of the con- 
ventional 2-flavor color superconductor, the spectrum of the fermionic quasi- 
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particles is different. In particular, two out of four gaped quasiparticles of the 
conventional 2SC phase become gapless in the g2SC phase. In addition, the 
number densities of the pairing quarks in g2SC at zero temperature are not 
equal [12]. For example, the density of red (green) up quarks is different from 
the density of green (red) down quarks. This is in contrast to the conventional 
"enforced" pairing [19] where the corresponding densities are equal. 

It is possible that the g2SC quark matter exists inside central regions of com- 
pact stars. This would be quite natural since this phase of matter is neutral 
with respect to electrical and color charges, and satisfies the /3-equilibrium 
condition by construction. Also, this phase is more stable than normal quark 
matter under similar conditions [12]. While suggestive, these arguments alone 
are not sufficient to make a conclusive statement regarding the chances of 
the g2SC phase to appear inside stars. In order to study this issue properly, 
one should consider the competition of the g2SC phase with all other possi- 
ble phases of quark/hadronic matter. However, the limited knowledge of the 
strongly coupled dynamics in QCD at realistic densities does not allow such 
a systematic quantitative study. 

In this paper, we study properties of the g2SC phase at finite temperature. 
Our interest in such properties is driven by the expectation that the g2SC 
phase of quark matter may already appear at a relatively early (protoneutron) 
stage of the compact star evolution when the temperature of the stellar core 
is quite high. If this is the case, the knowledge of finite temperature properties 
of quark matter would be needed to understand the stellar properties and, 
possibly, even to predict observational signatures of the new state of matter 
inside stars. 

This paper is organized as follows. In the next section, we set our notation 
by presenting a model of quark matter, and spell out the main assumptions 
in the analysis. Then, in Sec. 3 we discuss in detail the solution to the gap 
equation and the condition of the charge neutrality in the two-flavor quark 
model at zero temperature. In addition to reproducing the numerical results 
of our original short Letter [12], we also present many analytical results that 
have not been included in Ref . [12] . At zero temperature, the screening prop- 
erties of gluons in the g2SC phase, as well as the possible appearance of two 
additional spin-1 condensates, are deduced from the qualitative nature of the 
low energy spectrum of quasiparticles. The finite temperature properties of 
the g2SC phase are studied in Sec. 4. Among main results, we find that (i) 
the phase transition between the g2SC phase and normal quark matter is a 
second order phase transition, (ii) the ratio of the critical temperature and 
the zero temperature gap is not a universal number in the gapless phase, and 
(iii) a nonzero finite temperature gap may appear even if its zero temperature 
value is vanishing. Sec. 5 contains our conclusions and a brief discussion of the 
properties of the g2SC phase as well as their applications. 
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2 Model 



2. 1 Lagrangian density and parameters of the model 



Let us start our analysis by discussing the quark model that we use in the 
following. We assume that the strange quark is sufficiently heavy and does not 
appear at intermediate baryon densities under consideration (e.g., this might 
correspond to the quark chemical potential // = Hb/^ in the range between 
about 350 and 450 MeV). Therefore, we use the simplest SU(2) Nambu-Jona- 
Lasinio (NJL) model that was proposed in Ref . [20] . The explicit form of the 
Lagrangian density reads: 



L^q{iYd^ - mo)q + Gg [{qqf + {qil-,rqf 

+Gd [{ifee''-f^q){iqee%q'')\ , (1) 

where q^ = C(f is the charge-conjugate spinor and C = 27^7° is the charge 
conjugation matrix. The quark field q = qia is a four-component Dirac spinor 
that carries flavor {i — 1,2) and color (a — 1,2,3) indices. The Pauh matrices 
are denoted by r = (T\r^,r^), while (e)*'' = e^'^ and (e^)°^ = e'^^'^ are the 
antisymmetric tensors in the flavor and color spaces, respectively. We also 
introduce a momentum cutoff A, and two independent coupling constants in 
the scalar quark-antiquark and scalar diquark channels, Gs and Gd- 

The values of the parameters in the NJL model are chosen by fltting the pion 
decay constant /^r = 93 MeV and the chiral condensate (uu)^^'^ = (dd)^^^ = 
—250 MeV. In the case of mo 7^ 0, the additional parameter could also be flxed 
by fltting the value of the pion mass. In the rest of this paper, however, we 
consider only the chiral limit with toq = 0. Then the two model parameters, 
i.e., the coupling constant and the cutoff, are chosen as follows: Gs — 5.0163 
GeV"^ and A = 653.3 MeV [20,21]. Without loss of generality, we set the 
strength of the diquark coupling to be proportional to the quark-antiquark 
coupling constant, i.e., Gd = TjGs with a typical number for r] being around 
0.75 [12]. 

When considering the chiral limit in the following, we shall find that the con- 
stituent quark mass is zero in the model at hand in the g2SC phase. In reality, 
however, this may not be the case. At the same time, since the appearance of 
the g2SC phase is a Fermi surface phenomenon, our assumption should not be 
a strong limitation. Of course, nonzero quark masses should affect the Fermi 
momenta of different quark species. In turn, this could change slightly the 
mismatch parameter in our model (see below), and introduce some minor 
changes in the analysis [22]. However, to large extent, the results would not 
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change. 



2.2 Thermodynamic potential in (3 -equilibrium 



In /9-equilibrium, the diagonal matrix of quark chemical potentials is given in 
terms of quark, electrical and color chemical potentials, 

2 

/^y,a/3 = — IJ'eQij)Sal3 + -^y^s5ij{Ts)aP, (2) 



where Q and Tg are generators of U(l)em of electromagnetism and the U(l)8 
subgroup of the color gauge group. The explicit expressions for the quark 
chemical potentials read 



l^ur — l^ug — A* ~ gA*e + gA*8, (3) 
11 

IJ-dr = l^dg = g/^e + (4) 

2 2 

f^Ub^ f^- g/^e - (5) 

l^db= g/^e - (6) 

One should notice that, in general, there are two mutually commuting color 
charges related to the generators T3 and Tg of the SU(3)c group. Therefore, one 
could introduce two chemical potentials for these two different color charges. 
However, we require that quark matter in the 2SC ground state remains in- 
variant under the SU(2)c color gauge subgroup. This condition makes the 
introduction of the second nontrivial color chemical potential ix^ unnecessary. 

In the mean field approximation, the finite temperature effective potential for 
quark matter in /3-equilibrium with electrons takes the form: 



+ 



AG 



D 



127r2 

-E 



15 



AG, 



(2^ 



£;„ + 2Tln (l + e-^°/^) 



(7) 



where /i^ is the electron chemical potential and Qq is a constant added to make 
the pressure of the vacuum zero. For simplicity, we substituted the zero value 
of the electron mass which is sufficient for the purposes of the current study. 
The sum in the second line of Eq. (7) runs over all (6 quark and 6 antiquark) 
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quasiparticles. The explicit dispersion relations and the degeneracy factors of 
the quasiparticles read 



^± = ± [xl] (8) 

E^,^E{p)±fi,,, [xl] (9) 

E^^^E^ip)±Sfi. [x2] (10) 

Here we introduced the following shorthand notation: 



E(p) = ^p2 + m2, (11) 
E^{p)^J[E{p)±jl]^ + A\ (12) 



2 ~ 2 ~^~~6 Y' ^ ^ 

2 ~ 2 ~ Y' ^ ^ 

The knowledge of the thermodynamic potential is sufficient to calculate all 
other thermodynamic properties of quark matter in equilibrium. In the model 
at hand, the physical thermodynamic potential that determines the pressure of 
quark matter, Qphys = —P, is obtained from Q in Eq. (7) after substituting ^s, 
lie, m and A that solve the color and electrical charge neutrality conditions, 
i.e., 

on dQ , , 

ns = = 0, and uq = -—— = 0, (15) 



as well as the gap equations, 

_ = 0, and ^ = 0, (16) 



respectively. In the following two subsections, we shall derive the explicit ex- 
pressions for the charge neutrality conditions and the gap equations. 



2.3 Electrical and color charge densities 



By enforcing only the charge neutrality conditions in Eq. (15), one can con- 
struct the effective potential of the neutral quark matter as a function of two 
order parameters m and A. The subsequent minimization of such a potential 
determines the physical values of the Dirac mass of quarks and the value of 
the superconducting gap in the ground state of neutral quark matter. 
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From the potential in Eq. (7), we derive the expressions for the electrical 
charge density nq and the color charge density by taking partial derivatives 
with respect to and Hs, respectively. As a result, we arrive at the following 
explicit expressions: 



^8 = 



4 r (fp 



E-Ji 
EX 



El 



[l-/>A+)-/>A-)] 



+f{E:,) - f{E-,) + f{El) - f{Ej,) 



(17) 



and 



(27r)3 

/(^A+) - fiK^) + fiK-) + hEi-) 



(18) 



Here the standard Fermi distribution function has been introduced, f{E) ~ 
1/[1 + exp(£'/T)]. From Eq. (18) we notice that the hnear combination of 
charge densities hq + ^ng is given by a rather simple expression. In this partic- 
ular combination, the electrical and color charges of blue down quarks cancel. 
Had the dispersion relations in Eq. (10) been dependent on p, only, the charge 
cancellation from the corresponding quasiparticles would also be exact. It is 
the presence of the mismatch parameter S/j, that leads to a nonvanishing con- 
tribution. As we shall see, the latter survives even at zero temperature in the 
g2SC phase {S/i > A), while disappear in the 2SC phase {S/i < A). 

In the following, we shall require both charge densities in Eqs. (17) and (18) to 
be zero. Such a requirement is necessary in order to have a realistic construc- 
tion of dense quark matter that could appear in the bulk of a compact star. 
Indeed, without enforcing the neutrality, the repulsive Coulomb force would 
clearly overwhelm the gravity preventing the formation of a stable star. 

Before concluding this subsection, we mention that it also makes sense to 
study quark matter at fixed values of chemical potentials without imposing the 
conditions of local neutrality. Such studies are particularly appropriate when 
globally neutral mixed phases are considered [16,17]. In this paper, however, 
we do not discuss mixed phases, but concentrate our attention exclusively on 
locally neutral (homogeneous) phases of quark matter in /3-equilibrium. 
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2-4 Gap equations for order parameters 



Starting from the potential in Eq. (7), we derive an explicit form of the gap 
equations. In particular, the equation for the Dirac mass reads 



m — tuq = AGs 



I 



d^p m 



2^[l-/>A+)-/(^Z-)] 



+2^[i-/>i+)-/>i-)] 



(19) 



In the chiral limit (mo = 0), the trivial solution m — corresponds to a 
chirally symmetric phase of quark matter, while a nontrivial solution m ^ 
corresponds to a phase with spontaneously broken chiral symmetry. 

Similarly, we derive the gap equation for the diquark condensate. 



A = 8G' 



f d^p 
"^J (2^ 



A 



E 



^[l-/(^X+)-/(^A-)] 



A 



El 



[1 - HE, 



(20) 



This equation can also have trivial as well as nontrivial solutions for A. In the 
latter case, the corresponding ground state is a color superconductor. 



2. 5 Number densities of quarks 



According to one of the criteria of the g2SC phase, the densities of the quark 
species that participate in pairing dynamics are not equal at zero temperature 
[12]. This is in contrast to regular pairing in the conventional gaped color 
superconductors [19]. 

Even before studying the effect of the Cooper pairing on the quark densities, 
we should discuss how to derive the formal expressions for these densities. One 
appreciates immediately, of course, that the definitions of such densities are 
unambiguous only if the chemical potentials of all six types (2 fiavors times 3 
colors) of quarks are treated as independent parameters. This is exactly how 
we proceed. The quark densities are obtained by taking the partial derivatives 
of the potential fl in Eq. (7) with respect to these independent parameters. 
At the end, the relations between the chemical potentials required by the 
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/3-equilibrium arc imposed. The only point that needs some additional expla- 
nation is how to interpret the dependence of the thermodynamic potential in 
Eq. (7) on the parameters jl and Sjj,, entering through the quasiparticle ener- 
gies in Eq. (10). It turns out that p, and S/i should be replaced by ^{fJ^dg + f^ur) 
and \{lJ'dg — IJ-ur) 1 respectively, in one set of the double degenerate modes, and 
by \{^idr + A*uff) and \{^idr — A*ug) in the other set, see Eqs. (13) and (14). 

Because of the residual SU (2)c symmetry in the ground state of quark matter, 
the densities of the same flavor, red and green quarks are equal in the g2SC 
ground state. For example, the densities of the up quarks participating in the 
Cooper pairing read 



(27r)3 
E-p 
EX 



^[^-hEl,)-f{El.)] 

[l-/>X+)-/(^A-)] 



(21) 



Similarly, for the densities of the down quarks, we get 



(27r)3 
E-p 
EX 



:i-/>A+)-/(^A-)] 



-/(^A+) - ReU) + f{El-) + ~f{El_) 



(22) 



As we see, the densities of the up and down quarks participating in the Cooper 
pairing are not equal. In fact, the difference of the densities is given by 



^dg '^ur — ^dr '^ug 

(2^ 



/(E^,)+/(E+_)-/(^A+; 



(23) 



which is always nonzero at finite temperature provided the mismatch param- 
eter Sfi is nonzero. It is even more important for us here that this difference 
is nonzero at zero temperature in the g2SC phase of quark matter. 



n 



=0 iP'dr ^ 



ug)\T=0 



e {5ix - A) - + {5lif - A^) . (24) 
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This is in contrast to the 2SC phase {Sfi < A) where this difference is zero in 
agreement with the arguments of Ref. [19]. 

The unpaired blue quarks, Ub and df,, are singlet states with respect to SU(2)c 
symmetry of the ground state. The densities of these quarks are 



riub 



d' 



P 



(27r)3 

for the up quarks, and 



3 



(25) 



ridb = 2 



d^p 



(27r)3 

for the down quarks, respectively. 



3 I TT^ , 



(26) 



3 Gapless 2SC at zero temperature 



Here we present a detailed study of the zero temperature properties of dense 
two- flavor quark matter. A brief outline of this study was presented in Ref. [12]. 
In this paper, we include additionally many new results (e.g., regarding the 
gluon properties in g2SC phase of quark matter, and the appearance of spin-1 
color superconducting gaps on top of the g2SC phase) that were outside the 
scope of the original publication [12]. 

At zero temperature, the effective potential (7) for quark matter in /3-equilibrium 
with electrons takes the following simple form [12]: 



Here, without loss of generality, we assume that mo = (chiral limit). In this 

case, the constituent quark mass is also zero in the color superconducting phase 
of matter at sufficiently large baryon chemical potential. Thus, we substitute 
m = below. After making use of the quasiparticle energies in Eqs. (8)- (10), 
we derive 

O-O /^e , l^tb 

A 2 7 

-2|^(^(p + /.)2 + A2 + ^(p-;,)2 + A2) 
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29 {5^1 - A) / ^{5^1 - ^{p-fiy + A^), 



(28) 



where /x^ = ± y (5/^)^ — A^. The appearance of the last term that contains 
the step function 9{6n — A) is remarkable. It can be traced back to the exis- 
tence of two degenerate gapless modes in the quasiparticle spectrum, whose 
energies change sign at the values of momenta p — jj,^, see Eq. (10). In a 
different language, this sign change modifies the prescription for treating the 
quasiparticle poles in the complex energy plane. This is exactly the reason 
why the arguments of Ref. [19] regarding the conventional "enforced" pair- 
ing in color superconductors do not apply in the case of g2SC. The nonequal 
densities of pairing quarks, see Eq. (24), is the formal consequence of this fact. 



3.1 Solving neutrality conditions 

Now, let us proceed to the analysis of the color and electrical charge neutrality 
conditions in the model of dense quark matter at hand. At zero temperature, 
the corresponding two independent conditions could be written in the form: 
ns — and uq + ^Us = 0, i.e.. 



where the ellipsis denotes the terms of higher order in powers of A, and 




) 



) 



(29) 



riQ + 





+ {Sf, - A) ^{Sf^r - A2 + (M' - A^) 



= 0. 



(30) 
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A (MeV) 




(MeV) 



Fig. 1. The graphical representation of the solution to the charge neutrality con- 
ditions (thick dash-dotted line) and the solution to the gap equation for three dif- 
ferent values of the diquark coupling constant (thick solid and dashed lines). The 
intersection points represent the solutions to both. The thin solid line divides two 
qualitatively different regions, A < (5// and A > Sii. The results are plotted for 
= 400 MeV and three values of diquark coupling constant Gd = ^Gs with 
r] = 0.5, T] = 0.75, and r] = 1.0. 

This last condition is particularly simple in the 2SC phase, A > S/j,. In this 
case, the exact solution to this equation is available: 

3/i-2/i8 ^ for A>5fx. (31) 



The other equation, then, has the following approximate solution: 



45 /4 , 4^/3 \ 
^8^44 (5-6'/'-^)/^ + ---^ 0.0854/i+---, (32) 

where the ellipsis again denotes the terms of higher order in powers of A//i. 
As we see, the color chemical potential /ig is very small compared to the other 
chemical potentials, fie and /i. Of course, the corrections to the above solutions 
may become large with increasing the value of the gap A. Therefore, strictly 
speaking, these formulas are reliable only in the weakly coupled regime when 
A ^ /i. Our numerical results show, however, that /xg remains also small at 
relatively large values of the gap. Moreover, the approximation /zg works 
well in the region A < 6fi and in a wide range of parameters in the model. In 
the following, we shall use this observation to derive many useful analytical 
results. 
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Now, let us consider the neutrality condition in the region A < Sjj,. In our 
analytical treatment, we substitute /i^ = 0. As we mentioned above, this is a 
very good approximation. In the case of small values of A this should not be 
surprizing since the color charge vanishes in the normal phase when iig — 0. 
The value of the electrical chemical potential is determined by solving the 
neutrality condition in Eq. (30). The corresponding result for A) 
could be given by the inverse representation of the gap in terms of the potential 



^(nQ=0) 



This analytical result is barely distinguishable from the numerical solution rep- 
resented by the dash-dotted line in Fig. 1. Quark matter is positively charged 
to the left from the neutrality line, and it is negatively charged to the right. 

From Eq. (33), we find an approximate point of the intersection of the neutral- 
ity hne with the //g-axis (A = 0). This is denoted by a triangle in Fig. 1. This 
point determines the value of the electrical chemical potential in the normal 
phase of neutral quark matter. Its approximate value is given by a solution to 
a cubic equation, obtained from Eq. (33) after substituting A = on the left 
hand side. This value is fi'^^ ~ 0.2196/i. For example, at = 400 MeV, it is 
87.84 MeV which is in good agreement with the numerical result plotted in 
Fig. 1. 

Also, from Eq. (33), we determine the intersection of the neutrality line with 
the line A = dfi that separates two qualitatively different regions in the plane 
(/ie,A). This occurs at ~ |/i that corresponds to A*^^) = ^f^P — u)/^- 
For the choice of the quark chemical potential jj, = 400 MeV, used in Fig. 1, 
the corresponding values are //^^^ ~ 240 MeV and A^^^ ~ 120 MeV which are 
again in good agreement with the numerical result (cf., the location of the 
"knee" of the dash-dotted line in Fig. 1, denoted by a star). 

Before concluding this subsection, we should emphasize that the charge den- 
sities in Eqs. (29) and (30) do not depend explicitly on the diquark coupling 
constant. This is the property of the mean field approximation that would 
break after taking higher order corrections into account. It is reasonable to 
assume that the corresponding dependence is weak even in the full theory. As 
we shall see below, this property is not shared by the solution for the gap A 
which is very sensitive to changes in the coupling constant. This observation 
would have far reaching consequences below. 
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3.2 Solving gap equation 



Now, we turn to the solution of the gap equation (20) at zero temperature. In 
this case, the equation simphfies considerably. Its explicit form reads 



+ (^/X 2^ j m ^2 + 2(A2 - //2)2 + 

-e {5ii - A) [ 5iisJ{5iif - A2 + - i A^) In 



■ (34) 



As it should be, this equation always has the trivial solution A = 0. This 
corresponds to the normal phase of quark matter. If there is also a nontrivial 
solution A 7^ 0, the normal phase is not necessarily the ground state. In fact, 
the true ground state corresponds to the solution that minimizes the value of 
the thermodynamic potential (i.e., it maximizes the pressure). Moreover, in 
the case of an additional requirement of neutrality, the minimization should 
be done along the neutrality line discussed in the preceding subsection. 

For a moment, let us drop the condition of the charge neutrality and consider 
all possible nontrivial solutions to the gap equation (34). Wc realize right away 
that there are two different branches of solutions in two regions, A < S/i and 
A > S/j,, see Fig. 1. 

First, we start with the region A > 6^. This is the case of a relatively small 
mismatch between the Fermi momenta of the pairing quarks. In particular, 
it describes the conventional 2SC phase of quark matter. From Eq. (34), we 
derive the approximate analytical expression for the nontrivial solution, 

A~^(/x), for A>5/x, (35) 



where 



«,).2;X^e.p(-^ + ^). (36) 



Since the average chemical potential p, depends on /ig, see Eq. (13), the above 
solution for the gap is also a function of fi^.- This dependence, however, is 
very weak because /Xg enters /i with a relatively small weight. This is also 
confirmed by our numerical results for three different values of the diquark 
coupling constant presented in Fig. 1. With increasing /Xg, the value of the gap 
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decreases by less than 10% along the whole upper branch of the solution. In 
contrast, the dependence of the gap on the coupling constant is very sensitive. 

In the region A < S/i, the mismatch parameter S/i is large. This is where the 
g2SC phase of quark matter appears. The solution to the gap equation in this 
region is approximated well by the following expression: 



where ^ is the same weakly dependent function of //e that was defined in 
Eq. (36). By making use of this analytical expression, we see that the lower 
branch of the solution to the gap equation extends only over a finite range of 
chemical potentials, yU^™™) < Aig < fi^^^^^ (see Fig. 1). The values of fi^^"^^ and 
^(max) g^j^g determined by the solutions to the equations /ig = ^ifJ'e) and fie = 
2^(/ie). The corresponding graphical solutions are represented by the square 
and the diamond in Fig. 1, respectively. At fj, — 400 MeV and r] = 0.75, the 
corresponding approximate values are a; 107.9 MeV and /x^™'^'') pa 205.3 

MeV which are in reasonable agreement with numerical results in Fig. 1. 

In summary of this subsection, we note that the general nontrivial solution to 
the gap equation consists of two different branches in two regions, A > S/i and 
A < Sfi. The upper is the main branch that extends down to /Xg = (i.e., the 
case without a mismatch in the Fermi surfaces of the up and down quarks). The 
lower branch appears only in a finite window of electrical chemical potentials. 
It merges with the upper one at a point on the line A = S/i. In Fig. 1, 
three solutions to the gap equation ai /j, — 400 MeV in three regimes with 
different coupling constants, corresponding to 77 = 0.5, rj — 0.75, and r] = 1.0, 
are shown. The three curves have qualitatively the same shape, but differ 
by overall scale factors. As we shall see below, this difference in the overall 
scales has an important physical consequence. In particular, at weak coupling 
iv ^ O-'''), the ground state of neutral quark matter is the normal phase, while 
neutral g2SC and 2SC phases are most stable at intermediate (0.7 ^ f] ^ 0.8) 
and strong {rj > 0.8) couplings, respectively. This is studied in detail in the 
following subsection. 



3.3 Ground state of neutral matter 

The ground state of neutral matter is determined by the location of the mini- 
mum of the effective potential as a function of the gap. It is understood that the 
minimization is performed after the neutrality condition is fulfilled. In differ- 
ent words, the consideration is restricted only to the neutrality (dash-dotted) 




for A < S/j,, 



(37) 
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line in Fig. 1. Mathematically, the necessary condition for the minimum reads 
[dA ~ dA diiAdA ^ dfis[dA ^ ^ 

\ / nQ,ns=0 f-e \ / nQ,ns=0 r~o \ / nQ,ns=0 



We notice that the last two terms with the derivatives taken along the neu- 
trality line are proportional to the electrical and color charge densities. They 
become zero when the requirement of neutrality is enforced. Therefore, the 
above condition formally coincides with the gap equation studied in the pre- 
ceding subsection. The important difference is that Eq. (38) is defined only at 
the neutrality line in the {/ig, A) plane. 

In order to satisfy the condition in Eq. (38), it is very convenient to use a 
simple graphical construction, presented in Fig. 1. The extremum values of 
the gap are determined by the points of intersection of the neutrality (dash- 
dotted) line with the solutions to the gap equation (i.e., either with one of the 
branches of the non-trivial solution, or with A = line). 

At very weak coupling (e.g., rj = 0.5 case), both branches of the nontrivial 
solutions to the gap equation are squeezed into a small region around the 
origin in the (/ig, A) plane. As a result, the neutrality line does not intersect 
them. Therefore, the only extremum of the effective potential corresponds to 
the trivial solution A = 0. As is easy to see, this is the global minimum. Thus, 
the normal phase is the ground state of quark matter in this regime. 

In the opposite case of very strong coupling (e.g., rj — 1 case), the potential 
has two extrema. Of course, one of them corresponds to the normal phase 
(A = 0). This is a local maximum of the potential. The other extremum is a 
minimum that corresponds to the 2SC phase with A > 6fi. Note that the types 
of the extrema could be inferred from the expected topology of the effective 
potential after taking into account that this potential is bounded below when 
A ^ oo. Our direct calculation show that the 2SC phase is indeed the ground 
state. At large coupling, this is hardly a surprise. 

The regime of intermediate coupling (see, e.g., rj — 0.75 case) is most inter- 
esting. At first sight, it looks similar to the case of strong coupling because 
the potential has two extrema. One of them occurs at A = and the other at 
A ^ 0. The peculiar thing is that the nontrivial extremum corresponds to an 
intersection of the neutrality line with the lower branch (defined in the region 
A < 5/x) of the nontrivial solutions to the gap equation. By making use of 
the same topological arguments as in the previous case, we conclude that the 
nontrivial solution defines the ground state of neutral quark matter. This can 
also be checked by a direct numerical calculation of the effective potential, see 
the solid line in Fig. 2. As we shall discuss later in more detail, this ground 
state with A 7^ is the g2SC, rather than ordinary 2SC phase. 
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Fig. 2. The effective potential as a function of the diquark gap A calculated at a 

fixed value of the electrical chemical potential /ig = 148.4 MeV (dashed line), and 
the effective potential defined along the neutrality line (solid line). The results are 
plotted for fx = 400 MeV and Go = vGs with r) = 0.75. 

At this point, it is relevant to make a few comments regarding the early studies 
of the color superconducting quark matter at finite values of the electrical (or, 
equivalently, isospin) chemical potential [23]. There, the electrical chemical 
potential was treated as a free parameter, and the charge neutrality condition 
was not imposed. Obviously, the 2SC phase is the ground state of quark matter 
at small values of the mismatch parameter S/i (which is identical with 
here). However, when the value oi S/j, increases and reaches a critical value, 
a first order phase transition happens. The approximate value of the critical 
electrical chemical potential is fif ^ V^Ao where Aq is the value of the gap 
at /ig = (strictly speaking, this estimate is derived in the approximation 
of weak coupling [24]). At higher values of 6n, the normal phase is the new 
ground state of quark matter. In this picture, there is no room for the g2SC 
phase. 

Now, how can one reconcile the appearance of the g2SC phase after imposing 
the charge neutrality conditions, and its absence without such conditions? To 
answer this, we refer again to Fig. 1. In absence of the neutrality condition, one 
defines the potential V^( A) by restricting J7(A, fi^) to a vertical line He — Const 
for each fixed value of the parameter /Xg. As a result, one ends up with a whole 
family of the effective potentials that are functions of A. By changing the 
parameter /ig, one goes from one representative of the family to another. The 
extrema of the effective potentials are determined by the intersection points 
of lie — Const line with the line of solutions to the gap equation. 
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Sweeping the values of ji^ from left to right in Fig. 1, one easily understands 
the reason for the first order phase transition in the quark matter. At small 
/^e, the potential has two extrema at A = and at A 7^ 0. Naturally, the latter 
is the ground state. At intermediate values of ji^ where the lower branch of 
the solution to the gap equation exists, //^™°^ < jie < A^^™^^ , the potential has 
already three extrema that are given by the intersection of the line He = Const 
with the /ig-axis (i.e., A = solution) and with the two branches of the 
nontrivial solutions. Interestingly, the solution on the lower branch corresponds 
to a maximum of the potential. The same solution in the neutral matter was 
the global minimum of the properly defined potential, see Fig. 2. The other 
two extrema are minima. It is the competition between them, that results in 
a first order transition from the 2SC to the normal phase. 

It should be clear that this first order phase transition is unphysical under the 
requirement of local neutrality of quark matter. This is because the transi- 
tion typically happens between two types of matter with different, in general, 
nonzero charge densities (e.g., positively charged color superconducting matter 
and negatively charged normal quark matter). It is worth mentioning, how- 
ever, that this first order phase transition can get physical meaning in globally 
neutral mixed phases of quark matter [16,17], and in some condensed matter 
systems where there is no analogue of the charge neutrality condition [25] . 



3.4 Weak, intermediate and strong coupling regimes 



Let us discuss the rigorous meaning of the weak, intermediate and strong 
coupling regimes in the model at hand (as we shall see below, this could get 
somewhat ambiguous at finite temperature). We start from the weak regime 
first. By definition, the normal phase is the only neutral phase in this case. As 
should be clear, this happens when there is no intersection of the neutrality 
line with the two nontrivial branches of the solution to the gap equation. 

From Eq. (33) and the discussion following it, we know that the neutrality 
line intersects the /ig-axis at fie = jJi^P ~ 0.22yU (the corresponding point is 
marked by a triangle in Fig. 1). The lower branch of the gap solution, on the 
other hand, starts at fi^ — ^, see Eq. (37) and the point marked by a square 
in Fig. 1. Therefore, these two lines do not intersect if ^ < iJtS}\ and only the 
normal phase is consistent with the condition of neutrality, see the rj = 0.5 
case in Fig. 1. Now, we remind that while /j,'^^ is essentially independent of 
the coupling constant, the value of ^, for the definition see Eq. (36), is very 
sensitive to the choice of the diquark coupling. Taking this into account, we 
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derive the following approximate condition of the weakly coupled regime: 



Gn<G 



(1) 

D 



4A2 - 12^K4/ifln 



, "weak", (39) 



with /^i ~ — ^/x^^-* ~ 0.96/i here. At = 400 MeV, the above estimate gives 
Tj < 0.684. 

At larger couphng, the ground state of neutral matter is a broken phase. 
However, there are still two qualitatively different regimes depending on the 
location of the intersection of the neutrality line with the gap solution. If this 
occurs in the region A < S/j,, the g2SC phase is the ground state (e.g., r] — 0.75 
case in Fig. 1). Otherwise, the ground state is 2SC (e.g., r] — 1 case in Fig. 1). 

The criterion for the strong coupling regime reads 2^ > /ij,^-* = 0.6/i, i.e., 
^ > 0.3/i. To derive this, we note that the value of the electrical chemical 
potential //g = 0.6// corresponds to the intersection of the neutrality line with 
the boundary A = (this is marked by a star in Fig. 1). If this value is 
smaller than = 2,^ at which the gap solution intersects the boundary line 
(this is marked by a diamond in Fig. 1), then the ground state is the 2SC 
phase. This is satisfied if the value of the coupling strength in the di quark 
channel is sufficiently strong, i.e.. 



Gd > G^S - 



TT 



4A2 - 12/ii + 4/iiln 



4(A2-/ii)/(/.(^V2)' 



"strong", (40) 



with ^2 — IJ- — |A*i^^ — 0.9/x here. At /i — 400 MeV, this estimate gives 
7] > 0.806. 

Finally, the case of intermediate coupling strength is reahzed when < 
Gd < Cg^ or alternatively, //^ ~ 0.22// < ^ ^ 0.3/( ~ ^^^2. In this case, 
the ground state is the g2SC phase. We check that, at fj, = 400 MeV, this 
is realized when 0.684 ^ r/ ^ 0.806. A representative solution in this regime 
(with 7] = 0.75) is shown in Fig. 1 by the thick solid line. 

We would hke to mention that the critical values of the couphng constant G^^^ 
and G^' in Eqs. (39) and (40) are functions of the quark chemical potential 
fi, although the corresponding dependence is weak. Nevertheless, the above 
qualitative picture remains the same independent of the actual value of the 
quark chemical potential. Namely, there are always two critical values of the 
coupling constant that separate the three qualitatively different regimes. 
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Fig. 3. The schematic modification of the red and green quark dispersion relations 
in the 2SC phase (left panel), and the resulting quasiparticle dispersion relations at 
low energies (right panel). 

3.5 Quasiparticle spectra in 2SC and g2SC phase 

In the previous subsections, we established that the stable ground state in 
neutral quark matter is the g2SC phase, provided the coupling constant is 
neither too weak nor too strong (G^^ < < G^d)- However, these calcu- 
lations alone may not be sufficient to understand the physical properties of 
such a phase. To compensate this hmitation, we find it very useful to study 
the nature of the fermionic quasiparticles in the g2SC phase. 

It is instructive to start with the case of the ordinary 2SC phase (A > 5/i). 
All four fermionic quasiparticles that result from the Cooper pairing of the 
red and green quarks are gaped (note that, with the conventional choice of 
the gap pointing in the anti-blue direction in color space, the spectra of the 
blue up and blue down quarks are not affected by the pairing dynamics). 

The corresponding spectra are shown schematically in Fig. 3. The left panel 
gives a qualitative picture of the dispersion relations of the pairing quarks, see 
Eq. (10). The dark (light) shaded regions below (above) E = line represent 
the filled (empty) quark states. By combining the lines lying above E' = 
with the refiections of the lines lying below = 0, one obtains the complete 
dispersion relations of the fermionic quasiparticles in the 2SC phase. The low 
energy slices of these relations are shown in the right panel by dashed lines. 
For comparison, the dispersion relations of quasiparticles in free theory are 
also shown in the same panel by solid lines. 

It is interesting to observe how the nature of quasiparticles changes gradually 
when one goes from the region of small momenta to the region of large mo- 
menta. From Fig. 3, we see that the quasiparticles with the larger (smaller) 
gap are nearly identical with the down (up) quark holes at small momenta, 
p -C /2, and resemble the up (down) quarks at large momenta, p^ fx. Around 
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Fig. 4. The schematic modification of the red and green quark dispersion relations 
in the g2SC phase (left panel), and the resulting quasiparticle dispersion relations 
at low energies (right panel). 



the middle point, p — fi, the quark spectrum is strongly modified by the ap- 
pearance of the gap. Thus, the quasiparticles represent presumably a mixture 
of different types of quarks and holes. 

Now, let us turn to the quasiparticle spectra in the g2SC phase. Graphically, 
these are shown in Fig. 4 which has the same layout as Fig. 3. The left panel 
shows the qualitative dispersion relations of quarks in the case of A < 6fi, see 
Eq. (10), while the right panel gives the low energy slice of the corresponding 
quasiparticle dispersion relations (dashed lines). 

As one would expect, far outside the pairing region, p p,, the quasiparti- 
cle dispersion relations are similar to those in the 2SC phase. Also, around 
p c::^ fi, the gaped (double-degenerate) quasiparticle resembles the (double- 
degenerate) quasiparticle with larger value of the gap in the 2SC phase. This 
is the quasiparticle that interpolates between the down quark holes at small 
momenta, p <^ fl, and the up quarks at large momenta, p ^ Ji. The other 
(double-degenerate) quasiparticle is gapless. The energy of this latter vanishes 



at two values of momenta p = and p = fi'^ where fi^ = fi± J — A^. 



The most remarkable property of the quasiparticle spectra in the g2SC phase, 
see the right panel in Fig. 4, is that the low energy excitations [E <^ 6fi — A) 
are very similar to those in the normal phase represented by solid lines. The 
only difference is that the values of the chemical potentials of the up and 
down quarks fiur = f^ug and fidg = fJ^dr ai'e replaced by the values fi~ and 
/i"*", respectively. This observation suggests, in particular, that the low energy 
(large distance scale) properties of the g2SC phase should look similar to 
those in the normal phase. In the following subsection, we use this argument 
to justify that the Meissner effect might be absent in the g2SC phase. 
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3.6 Gluon properties in g2SC phase 



The main purpose for studying the NJL model in Eq. (1) is to understand 
the properties of dense quark matter under conditions that may exist inside 
compact stars. Indeed, this model should be adequate for getting an insight 
into the pairing dynamics of quarks and into some thermodynamic properties 
of quark matter. However, the NJL model lacks gluons. As reflection of this, 
it possesses the global instead of gauged color symmetry. The consequence is 
that there appear five Nambu-Goldstone (NG) bosons in a ground state of the 
model when the color symmetry is broken. In QCD, there is no room for such 
NG bosons. 

Perhaps, the NJL model can be thought of as the low energy theory of QCD 
in which the gluons, as independent degrees of freedom, are integrated out. 
The gluons could be reintroduced back by gauging the color symmetry in 
the Lagrangian density in Eq. (1), providing a scmirigorous framework for 
studying the effect of the Cooper pairing on the physical properties of gluons. 

Gluon properties in the 2SC phase of dense quark matter without the mis- 
match of the Fermi momenta of the up and down quarks have been studied in 
detail [26,27]. A similar quantitative analysis could also be performed in the 
case of a nonzero mismatch parameter 5^. Technically, this is not easy, and it is 
left for a future study. Here we argue, nevertheless, that some basic properties 
of gluons could be understood even without such an explicit calculation. 

At large densities, gluons become quasiparticles with rather complicated dis- 
persion relations. Typically, the only propagating (i.e., weakly damped) modes 
of gluons in dense medium are plasmons (see, however, Rcf. [28] where other 
interesting modes of dense quark matter are discussed). As in the case of 
metals, the minimum value of the frequency of plasmons is called the plasma 
frequency. Since its value in dense quark matter is of order y^/x, the plas- 
mons could be viewed as "high energy" modes that exist not only in the normal 
phase, but also in the 2SC and g2SC phases. The reason is that the properties 
of these high energy modes should not change much by the modification of 
the low energy quasiparticle spectra in broken phases. 

It is clear that the gluon self-energy carries information not only about the 
propagating modes, but also about highly damped modes and various screen- 
ing effects. The Meissner effect and the Debye screening are among them. It 
is important to mention that these types of large distance scale properties are 
determined to large extent by the low energy excitations. Therefore, they are 
very sensitive to the (dis-) appearance of gaps in quasiparticle spectra. 

First, let us discuss the Debye screening which is the property of the gluons 
of the electric type (i.e., spatially longitudinal gluons). In the case of the 
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2SC phase, this is studied in Ref. [26]. It is shown that the three gluons of 
the unbroken SU{2)f.. unlike the other five gluons of the original SU{3)c, do 
not experience the Debye screening in the far infrared region. This is the 
consequence of the absence of gapless quasiparticles in the 2SC ground state 
that interact with the gluons of the unbroken SU (2)c. Now, this property is no 
longer valid in the gapless phase studied here. One should expect, therefore, 
that all gluons of the original SU{3)c are subject to the Debye screening. Of 
course, the values of the Debye masses for different gluons need not be the 
same. 

A similar line of arguments may also be given in relation to the Meissner 
effect that affects five magnetic gluons (i.e., spatially transverse gluons). In 
this case, it is the appearance of a nonzero gap that leads usually to Meissner 
masses for the gluons [26,27]. We remind that the three gluons of the unbroken 
SU{2)c do not experience the Meissner effect. It is reasonable to suggest that, 
even with a nonzero mismatch parameter, five magnetic gluons in the 2SC 
phase (A > 6fi) should still experience the Meissner effect. Indeed, all those 
quasiparticles that were gaped at = still remain gaped. 

In the case of the g2SC phase (A < there appear gapless quasiparticles 
in the low energy sector of the theory, see Fig. 4. If the value of the quantity 
6fi — A is not vanishingly small, the low energy excitations {E < 6fi — A) in 
the g2SC phase look similar to the excitations in the normal phase in which 
the values of the chemical potentials are shifted slightly. Moreover, even the 
values of the density of states at the corresponding "effective" Fermi surfaces, 
i.e., at p — fj,~ and p — fj,'^, remain similar. By making use of the dispersion 
relations in Eq. (10), we derive 



dN 
dE 



25/i(/i" 



p=ij.- 



dN 
dE 



2Sfj,{fj,~ 



(41) 



Taking this into account, we may conclude that the large distance properties 
of the magnetic gluons in the g2SC phase should be the same as in the normal 
phase. This suggests, in particular, that there is no color Meissner effect in such 
a phase. This would be very unusual because the diquark gap is nonzero, and 
the Anderson-Higgs mechanism is expected to occur. Sometimes, this latter 
is identified with the Meissner effect. We emphasize, however, that the two 
effects are not the same. Of course, the Anderson-Higgs mechanism is realized 
in the g2SC phase in the following sense. The five NG bosons, related to the 
breaking of the global color symmetry, are unphysical in the g2SC phase. In 
other words, they are "eaten" by the gluons. This conclusion can be reached 
by using the gauge freedom in QCD that allows to get rid of the NG bosons 
(the corresponding gauge is called unitary). 
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Here, it is also instructive to mention that, in contrast to models with unbroken 
Lorentz symmetry, the gauge bosons in dense (or hot) medium have propagat- 
ing (massive) longitudinal modes even in the normal phase. Thus, the intuitive 
picture, being popular in particle physics, that connects the Anderson-Higgs 
mechanism with the appearance of nonzero masses for gauge bosons is not 
directly generalized to dense (or hot) medium. In the case of the g2SC phase, 
it might not be adequate even when applied to the appearance of the Meissner 
mass for magnetic gluons. 



3.1 Spin-l condensation 

The fact that the low energy fermionic quasiparticles in the g2SC phase look 
nearly the same as in the normal phase has yet another consequence. As we 
discuss below, there should appear two additional spin-l condensates on top 
of the g2SC ground state. In its turn, this should lead to the opening of small 
gaps around the points p = ^~ and p = n'^. Thus, strictly speaking, there will 
be no gapless excitations in the quasiparticle spectrum at zero temperature. 

Our argument here is based on the following simple observation. Around the 
effective Fermi momenta p = ^~ and p = fi'^, the gapless quasiparticles are 
double degenerate. This is the reflection of the degeneracy with respect to the 
red and green colors. Then, by taking into account the high degeneracy of 
states at p = and p = fJ.^ and the attraction in the color antisymmetric 
channel, we conclude that the Cooper instabihty should develop in the g2SC 
ground state around both points p — /i" and p = This instability is 
removed spontaneously by the appearance of two "secondary" condensates 
of Cooper pairs. In view of the Pauli principle, these are spin-l condensates 
of the so-called polar phase. The spatial rotation symmetry is broken in the 
presence of such a condensation [8] . 

Since spin-l condensates are much weaker than the spin-0 condensates, this 
additional condensations add a fine structure on the background of the g2SC 
phase. In practical terms, it means that many thermodynamic quantities, such 
as the energy density, the pressure and the particle densities should change 
very little. At the same time, other properties, such as the specific heat, the 
magnetic field penetration depth and various transport coefficients could be 
modified considerably by the appearance of the spin-l condensates. For ex- 
ample, the presence of these condensates will result in the color as well as 
electromagnetic Meissner effect at sufficiently low temperatures [9] . 

Before concluding this subsection, we would like to mention that an additional 
spin-l condensation of the blue up and down quarks is not expected in the 
model at hand. Even if there is an attraction in a flavor antisymmetric channel 
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(e.g., provided by instantons), the Fermi momenta of the blue up and down 
quarks are badly mismatched. This is the general outcome of enforcing the 
/9-equilibrium condition and the local neutrality in nonstrange quark matter. 



4 Gapless 2SC at finite temperature 



In this section, we study the properties of the g2SC phase at finite temperature. 
As we shall see, the interplay of the neutrality condition and the solution to 
the gap equation produces many unusual properties of quark matter that have 
no analogue in well known systems. One of our most intriguing observation 
will be that the ratio of the critical temperature to the value of the gap at zero 
temperature is not a universal number in the g2SC phase. In fact, depending on 
the coupling constant, this ratio could be arbitrarily large. The other striking 
property is related to a nonmonotonic dependence of the order parameter on 
temperature. In the most extreme case, for example, a nonzero gap can appear 
at finite temperature even it was exactly zero at T = 0. 

In a superconducting system, when one increases the temperature at a given 
chemical potential, thermal motion will eventually break up the quark Cooper 
pairs. In the weakly interacting Bardeen-Copper-Schrieffer (BCS) theory, the 
transition between the superconducting and normal phases is usually of second 
order. The ratio of the critical temperature T^^^ to the zero temperature value 
of the gap Aq is a universal value [29] 



rBCS = ^ = V ^ ^^^^ 

where ~ 0.577 is the Euler constant. In the conventional 2SC phase of 
quark matter with equal densities of the up and down quarks, the ratio of the 
critical temperature to the zero temperature value of the gap is also the same 
as in the BCS theory [3]. In the spin-0 color fiavor locked phase as well as 
in the spin-1 color spin locked phase, on the other hand, this ratio is larger 
than BCS ratio by the factors 2^/^ and 2^/^, respectively. These deviations are 
related directly to the presence of two different types of quasiparticles with 
nonequal gaps [30]. 

This commonly accepted picture of the finite temperature effects in supercon- 
ducting phases changes drastically in the case of dense quark matter when 
the /5-equilibrium and the local neutrality conditions are enforced. Below, we 
study this in detail. 
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4-1 Gap equation and charge neutrality condition 

Here we consider the gap equation and charge neutrahty conditions in the 
g2SC phase at finite temperature. Because of additional technical comphca- 
tions appearing at finite temperature, most of the results that follow will be 
obtained by numerical computation. It is important to keep in mind, however, 
that our approach remains conceptually the same as that at zero temperature. 

Before studying the gap equation (20), we recall that the right hand side of 
this equation is a function of the three chemical potentials /i, /Je and /ig. The 
values of He and Hs are not the free parameters in neutral quark matter. They 
are determined by satisfying the two conditions in Eq. (15). In our analysis, on 
the other hand, it is very convenient to enforce only the color charge neutrahty 
by choosing the function //§(//, //g. A) properly. As for the electrical chemical 
potential He, it is treated as an independent parameter at intermediate stages 
of calculation. Only at the very end its value is adjusted to make the quark 
matter also electrically neutral. 

As in the case of zero temperature, the resulting color chemical potential 
/ig is much smaller than the other two chemical potentials in neutral quark 
matter. This, of course, is not surprizing. Small nonzero values of /ig [typically, 
/ig ~ A^//i] are required only because an induced color charge of the diquark 
condensate should be compensated. Therefore, the smallness of /ig is protected 
by the smallness of the order parameter. The numerical calculations show, in 
fact, that even the approximation /is — does not modify considerably the 
exact solutions to the gap equation and to the electrical charge neutrality 
condition. In practice, we either tabulate the function /ig(/i, /ig, A) for a given 
set of parameters, or determine it numerically in the vicinity of the ground 
state solution. 

The solutions to the gap equation (20) for several different values of tem- 
perature (T = 0, 20, 40, 50.2 and 60 MeV) are shown graphically in Fig. 5 
(solid lines). The values of the temperature are marked along the curves. The 
results are plotted in the (/Xe,A)-plane. In computation we kept the quark 
chemical potential fixed, /i = 400 MeV, and used the diquark coupling con- 
stant Gd = rjGs with rj — 0.75. This is the same choice that we used at zero 
temperature in Fig. 1. 

Now let us briefiy discuss the temperature dependence of the solution. It is 
not very surprizing that the shape of the graphical solution is smoothed with 

increasing the temperature. The same applies to the disappearance of the 
double-branch structure of the solution at a finite value of the temperature, 
Tg ~ 30 MeV. Indeed, in a model that treats the electrical chemical potential 
as a free parameter, this value of the temperature Tg marks the expected 
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Fig. 5. The solutions to the gap equation (solid lines) and the neutrality condition 
(dashed lines) calculated for several values of temperature. The results are plotted 
for n = 400 MeV and the diquark coupling constant Go = ilGs with ry = 0.75. In 
this case, « 50.2 MeV. 

switch of two regimes. Namely, while the phase transition controlled by the 
parameter (i.e., neutrality is not required) is a first order phase transition 
at T < Ts, it becomes a second order phase transition at T > T^. 

We also observe that the value of the electrical chemical potential at the 
point where the solution to the gap equation intersects with the /ig-axis has 
a nonmonotonic dependence on the temperature. With increasing tempera- 
ture, this value increases first and, after reaching some maximum value, starts 
to decrease and goes down to zero eventually. When the interplay with the 
neutrality condition is taken into account later, this simple property of the 
solution would produce rather unusual physical results. 

We note that the temperature dependence of the solutions to the gap equa- 
tion is very similar to the dependence found by Sarma in a solid state physics 
analogue of the g2SC phase [25]. Of course, there was no analogue of the neu- 
trality condition in the system studied in Ref. [25]. Therefore, the mentioned 
similarity does not extend to the complete analysis of the g2SC phase that 
follows. 

In Fig. 5, we also show the neutrality lines (dashed lines) for the same values 
of temperature (T = 0, 20, 40, 50.2 and 60 MeV). The convention is that the 
lengths of the dashes decreases with increasing the value of the temperature. 
As we see, with increasing temperature, the neutrality line gets steeper while 
its intersection point with the /ig-axis moves towards larger values of He- The 
points of intersection of these neutrality lines with the lines of solutions to the 
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Fig. 6. The thermodynamic potential of neutral quark matter as a function of the 
diquark gap A calculated for several values of temperature. 

gap equation, when they exist, are shown as welL 

As in the case of zero temperature, see Subsec. 3.3 and Fig. 2, we need to show 
that the points of intersections of the gap solutions with the corresponding 
neutrahty hues in Fig. 5 represent the ground state of the neutral quark mat- 
ter. To this end, we calculate the dependence of the thermodynamic potential 
on the value of the gap A. Since the charge neutrality condition is satisfied 
only along the neutrality line, we restrict the thermodynamic potential only 
to this line. The numerical results for several values of the temperature are 
shown in Fig. 6. 

The results at /i = 400 MeV and i] = 0.75 show that, for any T < ^ 50.2 
MeV, the thermodynamic potential has the global minimum away from the 
origin, meaning that the corresponding ground state develops a nonzero ex- 
pectation value of A(T). Needless to say that this is the same value that one 
extracts from the location of the geometrical construction in Fig. 5. The ex- 
pectation value A(T) disappears gradually when the temperature approaches 
Tc from below. This is an indication of the second order phase transition. At 
temperature higher than Tc, the ground state of neutral quark matter is the 
normal phase with A(T) = 0. 

4-2 Temperature dependence of the gap 

The temperature dependence of the gap is obtained by numerical solution of 
the gap equation (20) together with the two charge neutrality conditions in 
Eq. (15). Of course, this is equivalent to the geometrical construction used in 



28 



120 
100 



> 80 

0) 

S 

zL 60 

to 



40 
20 









- 


> 


""•«^.^A(T,^e=o) : 


— ^ ^ 




\_ 
\ 












1 \ 
1 \ 


g2SC 


2SC 


1 \ * 

1 g2SC \ i 



10 20 30 40 

T (MeV) 



50 



60 



Fig. 7. The temperature dependence of the diquark gap (solid line) and the value of 
6fi = ^e/2 (dashed line) in neutral quark matter. For comparison, the diquark gap 
in the model with /Xg = and = is also shown (dash-dotted line). The results 
are plotted for // = 400 MeV and r] = 0.75. 

Fig. 5, where the intersection points of two types of hnes determine the values 
of the gaps in the ground state. 



The typical results for the default choice of parameters /i = 400 MeV and 
rj = 0.75 are shown in Fig. 7. Both the values of the diquark gap (solid 
line) and the mismatch parameter 6fi = fj^ef^ (dashed line) are plotted. One 
very unusual property of the shown temperature dependence of the gap is the 
nonmonotonic behavior. Only at sufficiently high temperatures, the gap is a 
decreasing function. In the low temperature region, T < 10 MeV, however, it 
increases with temperature. For comparison, in the same figure, the diquark 
gap in the model with /ie = and /is = is also shown (dash-dotted line). 
This latter has the standard BCS shape. 



Another interesting thing regarding the temperature dependences in Fig. 7 
appears in the intermediate temperature region, 22.5 ^ T < 37 MeV. By 
comparing the values of A(T) and 6fi in this region, we see that the g2SC 
phase is replaced by a "transitional" 2SC phase there. Indeed, the energy 
spectrum of the quasiparticles even at finite temperature is determined by the 
same relations in Eqs. (8) and (10) that we used at zero temperature. When 
A > 6fi, the modes determined by Eq. (10) are gaped. Then, according to our 
standard classification, the ground state is the 2SC phase. 

It is fair to say, of course, that the qualitative difference of the g2SC and 2SC 
phases is not so striking at finite temperature as it is at zero temperature. This 
difference is particularly negligible in the region of interest where temperatures 
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Fig. 8. The temperature dependence of the diquark gap in neutral quark matter 
calculated for several values of the diquark coupling strength ry = Gd/Gs- 

22.5 ^ r < 37 MeV are considerably larger than the actual value of the smaller 
gap, A — By increasing the value of the coupling constant slightly, however, 
the transitional 2SC phase can be made much stronger and the window of 
intermediate temperatures can become considerably wider. In either case, we 
find it rather unusual that the g2SC phase of neutral quark matter is replaced 
by a transitional 2SC phase at intermediate temperatures which, at higher 
temperatures, is replaced by the g2SC phase again. 

It appears that the temperature dependence of the diquark gap is very sensitive 
to the choice of the diquark coupling strength rj = Gd/Gs in the model 
at hand. This is not surprising because the solution to the gap equation is 
very sensitive to this choice. The resulting interplay of the solution for A 
with the condition of charge neutrality, however, is very interesting. This is 
demonstrated by the plot of the temperature dependence of the diquark gap 
calculated for several values of the diquark coupling in Fig. 8. 

The most amazing are the results for weak coupling. It appears that the gap 
function could have sizable values at finite temperature even if it is exactly 
zero at zero temperature. This possibility comes about only because of the 
strong influence of the neutrality condition on the ground state preference in 
quark matter. Because of the thermal effects, the positive electrical charge of 
the diquark condensate is easier to accommodate at finite temperature. This 
opens a possibility of the Cooper pairing that is forbidden at zero temperature. 

We should mention that somewhat similar results for the temperature depen- 
dence of the gap were also obtained in Ref. [31] in a study of the asym- 
metric nuclear matter. The asymmetry parameter in nuclear matter a = 
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(Pn — Pp)/{Pn + Pp) (where p„ and pp are the ncTitron and proton densities, re- 
spectively) plays a role similar to the coupling strength rj in the quark matter 
system. 



4.3 Nonuniversal ratio Tc/ Ao 



In the preceding subsection, we saw that the temperature dependence of the 
gap in the g2SC phase is very different from the standard benchmark result 
in the BCS theory. One of the main properties of the BCS temperature de- 
pendence is a universal value of the ratio of the critical temperature T^. to 
the value of the gap at zero temperature Aq, see Eq. (42). It is instructive, 
therefore, to calculate the same quantity in the g2SC phase of neutral quark 
matter. 

Let us start from a simple exercise, and consider the results plotted in Fig. 7 
more carefully. First of all, we find that T^^^ ^ 62.06 MeV and Af^ = 109.4 
MeV when there is no mismatch parameters in the model (i.e., He — and 
/is = represented by dash-dotted hne). It is not very surprizing that the ratio 
of interest r2sc' ~ 0.567 is in agreement with the BCS result. Now, if we check 
the results for the gap function in the g2SC phase (solid line), we find that 
rpg2sc ^ 50 2 MeV and Af^^ ^ 68.2 MeV. Thus, the ratio is rg2sc ~ 0.7357 
which is considerably larger than the 2SC result. 

It appears that the real situation is even more interesting. The result for the 
ratio rg2sc ~ 0.7357 is not universal. In fact, its value depends very much on 
the diquark coupling constant and, moreover, it can even be arbitrarily large. 
This last statement may not be so unexpected if we recall the temperature 
dependences of the gap shown in Fig. 8. 

The numerical results for the ratio of the critical temperature to the zero tem- 
perature gap in the g2SC function of the diquark coupling strength 
1] = Gd/Gs arc plotted in Fig. 9. The dependence is shown for the most inter- 
esting range of values oirj = G^/Gs, 0.68 < ?7 < 0.81, which allows the g2SC 
stable ground state at zero temperature, see the discussion in Subsec. 3.4. 
When the coupling gets weaker in this range, the zero temperature gap van- 
ishes gradually. As we saw from Fig. 8, however, this does not mean that the 
critical temperature vanishes too. Therefore, the ratio of a finite value of Tc to 
the vanishing value of the gap can become arbitrarily large. In fact, it remains 
strictly infinite for a range of couplings. 
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Fig. 9. The ratio of the critical temperature to the zero temperature gap in neutral 
quark matter as a function of the coupling strength r] = Go/Gs- 

5 Conclusion 



In this paper, we presented a detailed study of the zero and finite temperature 
properties of the gapless color superconducting phase of neutral two-flavor 
quark matter in /3-equilibrium. The spectrum of quasiparticles in this phase 
has four gapless and only two gaped fermionic modes, see Eqs. (8)-(10) with 
A < /ie/2. This is in contrast to the conventional 2SC phase that has only 
two gapless and four gaped quasiparticles, see the case with A > /ie/2. As we 
show, the choice of the ground state is predetermined by the strength of the 
diquark coupling. In strong coupling regime, the 2SC survives even under the 
constrain of the charge neutrality and /^-equilibrium. When the coupling is 
not very strong, the 2SC ground state is replaced by the g2SC phase. Finally, 
at weak coupling, the stable ground state corresponds to the normal phase. 

In the g2SC phase, like in the conventional 2SC phase, the original color 
symmetry SU(3)c is broken down to SU(2)c. Despite that, the low energy 
spectrum of the gapless phase looks very similar to the spectrum in the normal 
quark matter. By the "normal quark matter" here we understand the mixture 
of the following four normal components: (i) red and green up quarks with an 
effective Fermi momentum p^'^^ = fi~ , (ii) red and green down quarks with 
an effective Fermi momentum p'^'^^ = fi~^, (iii) blue up quarks with p"^ = fi^b 
and (iv) blue down quarks with pf^ = fidb- Needless to say that the similarity 
extends only in a small region of low energies, E < 6fi — A. This similarity 
suggests that there is no color Meissner effect in this phase of matter, although 
the gauge symmetry is broken through the Anderson-Higgs mechanism. 
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The same arguments regarding the resemblance of the low energy spectrum 
in the g2SC and normal phases of quark matter lead us to the conclusion 
that two "secondary" spin-1 condensates should develop on top of the gapless 
ground state. Indeed, the existence of highly degenerate gapless modes around 
the effective two Fermi surfaces should lead to the appearance of additional 
condensates, provided an attractive channel exists. Because of the red-green 
color degeneracy of the quasiparticles around p^^^ = ^~ and p^/"^ = the 
attractive antisymmetric channel should exist. As a result, two spin-1 gaps 
would open around the points where the gapless modes used to be. 

Most likely, the appearance of the spin-1 gaps is a small effect on top of the 
g2SC ground state. Therefore, many thermodynamical properties of quark 
matter would not be affected by their presence. The transport properties, on 
the other hand, should be modified. In particular, five gluons of the original 
SU(3)c gauge group should experience the color Meissner effect. In addition 
there will also be the electromagnetic Meissner effect [9] . 

At finite temperature, the properties of the g2SC phase are very unusual. One 
of the most striking results is that the ratio of the critical temperature to the 
gap at zero temperature is not a universal number. It depends on the choice 
of the coupling constant in the diquark channel. Moreover, the value of the 
ratio can become arbitrarily large with proper choice of the coupling. For a 
range of values, it is infinite. This is a simple consequence of the nonmonotonic 
temperature dependence of the diquark gap function on temperature. In an 
extreme case, the finite temperature gap can be nonzero even if the zero tem- 
perature gap is zero. To the best of our knowledge, the only other example of 
such a temperature behavior of an order parameter was obtained in Ref. [31] 
in a study of the asymmetric nuclear matter. 

In nature, gapless 2SC quark matter could exist in compact stars. Indeed, this 
phase is neutral with respect to electric and color charges and satisfies the 
/9-equilibrium condition by construction. If this is indeed the case, one may 
have a chance to detect the indirect signatures of its presence by deciphering 
the observational data from stars. The study of the physical properties of 
the g2SC phase in this paper is the first step in this direction. In future, one 
should also study the transport properties of the g2SC phase. For example, the 
knowledge of the neutrino emissivities and mean free paths would be crucial 
for the understanding of the cooling rates of stars with such a gapless phase 
in their interior. 

As we mentioned earlier, some stable gapless phases may also be realized in 
the asymmetric nuclear matter [31] and in two-component mixtures of cold 
fermionic atoms [15]. The results of this paper might also be useful for a deeper 
understanding of the properties of such nonrelativistic systems. For example, 
the similarity of the spectrum of the low energy quasiparticles in the gapless 
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phase and the normal matter seems to suggest that there is no superfluidity 
in the corresponding phases of cold atom mixtures. Of course, one should 
be careful before making a conclusive statement because the corresponding 
nonrelativistic systems are quite different from neutral quark matter. 
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